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Abstract. In this work we investigate the possibility to create cold Fi'2, RbFr, and 
CsFr molecules through photoassociation of cold atoms. Potential curves, permanent 

^ | and transition dipole moments for the Francium dimer and for the RbFr and RbCs 

molecules are determined for the first time. The Francium atom is modelled as a 
one valence electron moving in the field of the Fr + core, which is described by a 
new pseudopotential with averaged relativistic effects, and including effective core 

5_j ■ polarization potential. The molecular calculations are performed for both the ionic 

species Fr^~, RbFr + , CsFr + and the corresponding neutral, through the CIPSI quantum 
chemistry package where we used new extended gaussian basis sets for Rb, Cs, and 
Fr atoms. As no experimental data is available, we discuss our results by comparison 
with the Rb2, Cs2, and RbCs systems. The dipole moment of CsFr reveals an electron 
transfer yielding a Cs + Fr~ arrangement, while in all other mixed alkali pairs the 
electron is transferred towards the lighter species. Finally the perturbative treatment 
of the spin-orbit coupling at large distances predicts that in contrast with Rb2 and 
Cs2, no double- well excited potential should be present in Fr2, probably preventing an 
efficient formation of cold dimers via photoassociation of cold Francium atoms. 
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1. Introduction 

As discussed for instance by Sprouse et al pQ the recent development of laser cooling 
and trapping of radioactive atoms opens ways for new investigations like the search 
for the electric dipole moment (EDM), (3 decay, cold atom-atom collisions, Bose- 
Einstein condensation and more precise atomic clocks. As the heaviest alkali-metal 
atom, Francium offers a sensitive test of the influence of relativistic and quantum 
electro dynamic effects (QED) on the atomic structure. It is also an excellent system to 
study weak interactions in atoms such as atomic parity non-conservation (PNC) J2HH]- 

However, none of the Francium istopes is stable. The longest-lived isotope is 223 Fr 
with a half-time of 21.8 min, while the shortest-lived isotope 215 Fr has a half-time of 0.12 
/is. Several groups have been successful in trapping long-lived isotopes in a magneto- 
optical trap (MOT): in 1995, Simsarian et al captured 210 Fr atoms (half-time 3.2 min) 
produced at the Stony Brook LINAC |H E]; a collaboration between the groups of 
C. Wieman from the University of Colorado and H. Gould from Lawrence Berkeley 
Laboratory lead to the trapping of 221 Fr (half time 4.8 min) 6J; more recently, a 210 Fr 
MOT has been set up in Legnaro (Italy) [ZJ, but to our knoweledge its optimization is 
still under way. 

Even if the spectroscopy of Francium atom is not yet fully investigated, a 
noticeable amount of experimental and theoretical data is now available for energy 
levels, hyperfine splittings, and transition dipole moments. In contrast, experimental 
conditions allowing to form Francium dimers have not yet been achieved, and the 
availability of cold Francium atoms in MOT's may allow one day to perform the 
photoassociation spectroscopy of such radioactive molecules. The creation of a mixed 
species cold atom trap involving Francium and Rubidium or Cesium appears as a closer 
possibility [S], which would allow to explore the RbFr or CsFr heteronuclear molecules. 

On the theoretical side, no potential curves have been published yet on these 
molecular systems. To our knowledge, the only results concern the bond length and 
the dissociation energy of Fr2 which have been determined by Lim et al [Oj using various 
approaches based on relativistic coupled-cluster method or functional density theory. 
Also the van der Waals interactions between Fr, Rb, and Cs atoms have been computed 
by Derevianko et al |l()j . 

The purpose of this paper is to calculate the electronic properties of Fr 2 , RbFr, and 
CsFr and of their cations, including potential curves, permanent and transition dipole 
moments. Our aim is also to predict the rates for the photoassociation and for the 
formation of cold molecules, as we did previously for all mixed alkali dimers involving 
Li to Cs atoms [O]. Such predictions should be useful for upcoming experimental 
investigations. Our theoretical method relies on effective potentials for modelling the 
Fr + ionic core, as in our previous paper (hereafter identified as paper I) [12] devoted to 
the accurate determination of permanent dipole moments of heteronuclear alkali dimers. 

The combination of large-core pseudopotentials (namely one active valence electron 
per alkali atom), complemented by core-polarization pseudopotentials (CPP) following 
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the idea of Miiller and Meyer ^3] has proved up to now to be one of the most successful 
approaches to accurately deal with the electronic structure of alkali pairs, even up to 
cesium [HJ El EH El CHI UE] • One of the advantage of the large-core+CPP approach is 
that it allows full configuration interaction (CI) even in large basis sets for the valence 
electrons (reduced to 2 in the case of alkali diatomics), and that the core polarization 
mimics the effect of coupled core-valence double excitations that otherwise should be 
calculated in the CI. Moreover, even for compounds including cesium, most calculations 
were performed within the framework of the so-called scalar (i.e. spin-orbit averaged) 
pseudopotentials, or averaged relativistic pseudopotentials (AREP) [2011211122], which 
take into account scalar relativistic effects (mass correction, Darwin term) [23J. The 
definition of quasi-relativistic pseudopotentials is actually closely linked to the way 
electrostatic and spin-orbit coupling are treated. AREP's are very convenient for 
molecular calculations, as configuration interaction is achieved in a first step, whereas 
spin-orbit coupling responsible for multiplicity splitting is accounted for a-posteriori. 
Other schemes have been implemented in the case of heavy elements where spin-orbit 
coupling is large, namely two-component schemes [213 m ^ ne basis of spinors which treat 
electrostatic and spin-orbit coupling on an equal footing in the configuration interaction. 
Four-component CI calculations are also sometimes feasible J2H] but still require a heavy 
computational effort with respect to the other methods. Detailed and recent overviews 
of the various techniques can be found in references J2E1 123 EE1 123 EDI E01 EU IH2j- 
In alkali diatomics, spin-orbit coupling does not directly concern the ground state but 
only the excited states, and it therefore remains relatively small, even for lowest excited 
states of the heavy elements: the Cs(6p) fine structure splitting is 554 cm -1 , and the 
one of Fr(7p) is 1686 cm -1 [23], to be compared for instance with the one of iodine, 
close to 1 eV. Core-valence correlation is a larger effect in heavy alkalis, especially 
for the lowest valence states, closer to the core. Even more important is the two- 
electron valence correlation in alkali diatomics which must be treated accurately. This is 
why quasi-relativistic calculations with averaged relativistic large-core pseudopotentials 
complemented by CPP's proved to be very accurate for complexes of alkali atoms at 
least up to cesium [341 . 

Given this ability, it is appealing to extend this simple scheme to systems including 
the heaviest atom of the series, namely Francium. We have therefore determined a 
spin-orbit-averaged semi-local pseudopotential for francium in the large-core approach, 
namely with a single external electron, using the shape-consistent norm-conserving 
technique of Durand and Barthelat [3J3 [3H], complemented with core-polarization 
pseudopotentials (CPP). We describe in detail in section [2] the chosen model, that we 
compare to the model potential approach of Klapisch [37]. We recall in section [3J the 
main steps of the molecular calculations, and our results for molecular ions and neutral 
diatomics of Francium coumpounds. Atomic units will be used throughout the text, 
except otherwise stated. 
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2. Atomic structure calculations and pseudopotential determination 

In 1939, M. Perey of the Curie Institute in Paris discovered a new radioactive element 
that she named Francium [SH]. The optical spectrum of Francium remained unknown 
until 1978 when the D2 line of francium was detected by Liberman et al at CERN on the 
on-line mass separator ISOLDE JHH]. Then the precise spectroscopy of francium isotopes 
was performed by the same group still at CERN [HUES], while the properties of Fr highly 
excited states were explored at Moscow University |41j. After 1995, spectroscopic studies 
were carried on with cold trapped Francium atoms [H [HI E21 E] ■ Hyperfme splittings 
[HnjlI3| and radiative lifetime measurements [121 E] became then available. Among the 
energy levels, s levels and Rydberg d series have been the most investigated using two- 
step or two-photon excitation from the Ip manifold jUJ HUJ ESI Hi] , as well as the 7pj and 
8pj multiplets. An ionization potential of 32848. 972(9)cm~ 1 has been determined from 
the quantum defects of ns and nd Rydberg levels [33J. However the lowest 6d 3/2,5/2 levels 
have not been observed yet, even if their location has been predicted using quantum 
defects fits of ndj energies j46j . Levels of higher symmetry (/, g) as well as p levels 
beyond n = 7, 8 remain unknown. 

On the theoretical side, Dzuba et al jUlEHj used relativistic Hartree-Fock and many- 
body perturbation theory to calculate the transition energies, fine structure intervals, 
electric transition amplitudes, and addressed the issue of parity violation in this atom. 
Eliav et al jH| calculated ground and excited transition energies as well as ionization 
potentials and electron affinities [SHI IS] for all alkali atoms including Fr, with the 
relativistic coupled-cluster method. Relativistic many-body calculations of transition 
energies, hyperfine constants, electric-dipole matrix elements and static polarizabilites 
were performed by Safronava et al for all alkali atoms [H2]- Derevianko et al [SB] used 
many-body perturbation theory (MBPT) and model-potential calculations to investigate 
low-energy photoionization parameters for Fr. The same model-potential approach has 
been adopted by Marinescu [HI] to compute long-range parameters for all alkali pairs 
but Fr and by Marinescu et al jHHj for a Francium atom pair. More recently, long- 
range parameters for homonuclear and heteronuclear alkali dimers including Fr have 
been calculated by Derevianko et al jHH UHj using relativistic coupled-cluster method. 

All these data represent a rich starting point for our molecular calculations which 
first impliy the determination of an accurate effective potential for the representation 
of the Francium ionic core. We performed three different atomic calculations. We 
first implemented for Francium the parametric energy-adjusted model potential (MP) 
proposed by Klapisch [HZ|, further improved to involve ^-dependent parameters (with 
I the orbital moment of the valence electron), core polarization effects, and spin-orbit 
terms j27j. Accurate atomic energy levels are obtained, allowing the use of MP predicted 
levels to check the quality of the pseudopotential described below. However, the Klapisch 
type model potentials conserve the nodal structure of the wavefunctions even within 
the core region, which is not very convenient for molecular calculations, especially in 
such heavy elements. In line with previous works on alkali diatomics, we preferred 
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to use a pseudopotential approach which yields nodeless atomic wave functions in the 
core region. We thus adopted a mean-field approach extraction. We first performed 
relativistic Dirac-Fock SCF calculations (DF) with Desclaux's package [SB], to obtain 
single electron energies and spin-orbitals considered as reference data for subsequent 
pseudopotential determination. We then determined a spin-orbit-averaged semi-local 
pseudopotential (PP) for Francium within the large-core approach, complemented with 
core-polarization pseudopotential (CPP) which will be used in molecular calculations of 
section 01 

2.1. Parametric model potential (MP) 

Such a model is known to be efficient for instance to calculate multiphoton cross sections 
of heavy alkali atoms K and Cs |23 EH] and complicated spectra of heavy alkaline-earth 
atoms jSTj. A similar model potential without spin-orbit has been also used by other 
authors for the determination of the van der Waals interaction between alkali atoms [HI] , 
or of photoionization parameters in Fr [53] • The semi-local model potentiel accounting 
for the interaction between the ionic core and the external electron depends on the 
empirical parameters cq and contains a polarization term involving the static dipole 
polarizability (core- polarization contribution) of the ionic core a^: 

Vt(r) = -I{l+(Z-l)exp(-^r)W 2 rexp(-^r)}-|^{l-exp[-(r/rf) 6 ]}(l) 

where Z is the nuclear charge, and r e c are cut-off radii truncating the polarization 
potential at short electronic distances r. We use the theoretical value for ^^=20.38 au 
[o r 2] obtained with fully relativistic coupled cluster method, in very good agreement with 
the value 0^=20.41 au computed using relativistic random phase approximation [56J. 
The one-electron radial function P n £j (where j is the total angular momentum of the 
external electron) is the solution of the radial Schrodinger equation with eigenvalue e n ^\ 

r_l_^_ 1(1+1) 

^ 2dr^ + 2r 

The spin-orbit scalar operator V^{r) is explicitely included and is approximated 
by: 

where a is the fine structure constant. The last factor suggested by the Dirac 
equation [H3] is included, to ensure that the solutions of the radial Schrodinger equation 
are well-defined when r — ► 0. The empirical parameters a\ and r e c are adjusted until the 
energies e n gj agree with the experimental atomic energies (table H}. 

Atomic binding energies of Francium obtained by the MP approach are compared 
with experimental values and with several other computations in table 121 the MPBT 
calculations of Dzuba et al performed whithin the Bruckner expansion |48j . the 



a + "Vt^ + Vt(r) + V s f{r) - e n£j ]P nej = (2) 
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I 


a{ 


a{ 


<4 


i 

c 




1 

2 


4.29673059 
4.10987255 
3.39019850 


13.58990866 

5.87755911 

0.78543635 


1.54459081 
1.36939901 
0.92256617 


1.35534853 
1.22513294 
1.96472431 



Table 1. Empirical parameters involved in the model potential of eqQ] describing the 
Francium ionic core. The value ay =20.38 a.u. |62| is introduced in the polarization 
potential. 



Level 


Exp 


MP 
Present work 


AE 


48 


AE 


49 


AE 


£21 


AE 


|53|(1) 


E3|(2) 


7s 


-32848.87 


-32848.87 


0.05 


-32762 


86 


-32839 


9.8 


-32735 


113 






7p- 


-20611.46 


-20613.03 


-1.5 


-20654 


-42 


-20584 


27 


-20583 


28 






6d- 




-16684.56 




-16623 




-16370 












6d+ 




-16412.47 




-16423 




-16194 












7 P + 


-18924.87 


-18922.33 


2.5 


-18926 


-1.2 


-18913 


11 


-18907 


17 






8s 


-13116.62 


-13116.39 


0.2 


-13082 


35 


-13121 


-4 


-13051 


65 






8p- 


9735.91 


-9730.38 


5.55 


-9742 


-6 


-9719 


17 


-9712 


23 






8p+ 


9190.56.0 


-9199.62 


-9.0 


-9188 


-29 


-9158 


32 


-9176 


14 






7d- 


-8604.04 


-8605.25 


-1.2 


-8663 


-58 














7d+ 


-8515.57 


-8514.55 


1.0 


-8574 


-58 














9s 


-7177.85 


-7177.94 


-0.09 


-7160 


18 






-7148 


29 






9p- 




-5733.61 




-5736 








-5724 




-5748 


-5737 


9p+ 




-5493.70 




-5485 








-5480 




-5496 


-5487 


8d- 


-5248.22 


-5246.58 


1.6 


-5261 


-12 














8d+ 


-5203.50 


-5204.42 


-0.9 


-5218 


-15 














lOp- 




-3787.81 












-3782 




-3800 


-3790 


10p+ 




-3658.66 












-3650 




-3662 


-3655 


10s 


-4538.26 


-4538.29 


-0.03 


-4534 


4 






-4522 


16 







Table 2. Theoretical binding energies of Francium (in cm -1 ) obtained with the present 
MP approach, compared with other relativistic calculations and experimental energies 
from refs. |33l|3§]. Differences AE between calculated and experimental energies are 
also displayed. Levels labelled with a ± sign means j = £ ± 1/2. 



relativistic coupled-cluster (CC) calculations of Eliav et al 49J, the relativistic many- 
body calculations of Safronova et al |H2] , and the results of the Quantum Defect Theory 
(QDT) combined with MBPT and model potential calculations of Derevianko et al [53J. 
The present MP energies are found in good agreement with the available observed 
energies, and are generally more accurate than the other theoretical predictions. The 
fine structure splitting of the 6d level, AE so = 272.1 cm -1 , is found larger than in the CC 
value of ref.|49J, AE so = 200 cm -1 , and also disagrees with the values predicted using 
quantum defects fits of ndj energies jl^l {AE so = 188 cm -1 ). Our predictions for all 
the other unobserved npj levels agree with other theoretical values, giving in particular 
similar fine-structure splitting for 8p to lOp levels. These predictions for energy levels 
will be used in the next paragraph. 



Electronic properties of Francium diatomic compounds and prospects for cold molecule formation! 



1 


&i,e 


Ci,t 


rn,e 





0.784936 


14.194790 


1 


1 


0.170429 


0.352240 







0.170429 


0.309530 


1 


2 


0.258663 


-3.442104 







0.258663 


4.652569 


-1 



Table 3. Parameters of the pseudopotential of eq0] representing the Francium core, 
implemented in the present work. 



2.2. Pseudopotential approach (PP) 

We actually used this approach to perform the molecular calculations for Fr 2 and 
RbFr and CsFr compounds. The ionic core is described by a semi-local ^-dependent 
pseudopotentiel (PP) [H21 HIE] employed in paper I: 

W t (r) = -1/r + £ C hl r n ^ exp(-c^r 2 ) (4) 



As explained in paper I, and in contrast with the MP approach, the parameters Cm 
and an in equation|Ufor Rb and Cs ionic cores are adjusted to reproduce the energies and 
valence orbitals of an all-electron Hartree-Fock Self-Consistent Field (SCF) calculation 
for the atomic ground state. 

The Francium PP parameters (table 0} are chosen in order to match the results of 
the all-electron relativistic Dirac-Fock (DF) SCF calculations performed with Desclaux's 
package [SHI- More precisely, we determine the large components of the Dirac-Fock 
solutions for the valence levels 7si/ 2 , 7pi/ 2 ,3/2, 6^3/2,5/2 and 5/5/2,7/2, and we adjust the 
PP parameters to match the corresponding spin- orbit- averaged orbitals. Relativistic 
corrections to atomic radial wave functions such mass- velocity, and the Darwin terms 
are automatically introduced in Dirac-Fock orbitals and thus in the pseudopotential [23J, 
which keeps the simple expression of equationHJ As discussed by several authors on Cs or 
Fr exemples |4T | l47" | 152*] . zero-order Dirac-Fock energies do not include correlation effects 
necessary to reproduce the atomic energy level structure, taken in account through 
MPBT in ref. (13 0H "H2J EHj- However, while intra-core correlation effects are hardly 
accessible via CI, inter-shell core- valence correlation effects can be described via the core- 
polarization operator. Let us note finally that an ab initio small-core pseudopotential 
(involving an ionic core with 78 electrons, instead of 86 in the present work) including 
relativistic effects has been obtained by Ermler [64^ , and further used by Lee et al to 
compute ground state potential energy curves of MO oxydes with (M=Rb,Cs,Fr) [65J. 
However the authors have not extended their investigations towards francium-alkali 
pairs. 

The Dirac-Fock orbitals for the fine structure doublets 7pi/2,3/2, 6^3/2,5/2, and 
5/5/2,7/2 are plotted in figure "D It is clear that even if differences between 7pj orbitals 
are visible, this figure suggests that spin-orbit averaged orbitals can be used to set up 
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Figure 1. Dirac-Fock orbitals for the (a) 7pi/ 2 ,3/2, (b) 6d 3 / 2 ,5/2, (c) 5f 5 / 2j 7/2 levels 
of Francium. The orbital with the lowest (largest) j value is drawn with full (dotted) 
line. 

Figure 2. Spin averaged orbitals for the (a) 7s, (b) 7p, (c) 6d levels of Francium, 
calculated with various approaches: Dirac-Fock (dashed line), present pseudopotential 
(red full line), present model potential (blue dotted line). 

a pseudopotential for the Francium ionic core with averaged relativistic effects. Such 
orbitals are drawn in figure |2] for the 7s, 7p, and 6d levels, as computed via the present 
MP, DF, and PP approaches. As expected, the mean-field PP pseudo-orbitals perfectly 
match the DF spin-orbit averaged orbitals in the region outside the core, i.e. the outer 
lobe of the wave function is well reproduced. In contrast, the MP orbitals are shifted 
inwards with respect to the DF orbitals, except for the 7s wavefunction. The same 
trend (not shown here) is actually observed with the nodeless PP pseudo-orbitals when 
the core polarization effects (see next section for the modelling of these effects) are 
included, compared to those pseudo-orbitals without core polarization. As the MP 
orbitals intrinsically contain core polarization effects, we suspect that this explains their 
observed pattern. 

2. 3. PP results for atoms 

The external electron in each atom (Rb, Cs, or Fr) is described with a Gaussian basis set. 
In order to check the accuracy against the size of the basis set, we proceed as in paper 
I, defining two basis sets for Rb and Cs: we first use a [7s4p5<i] contracted to a [6s4p4<i] 
basis, identical to the one introduced by Pavolini et al jHj previously refered to as basis 
"A" in paper I. Then we consider for Rb the large [9s6p6d4f] basis of uncontracted 
Gaussian functions resulting in 138 molecular orbitals (basis "B" in paper I), and we 
elaborate a similar [9s6p6d4:f] basis for Cs refered to as "B"' in the following (table EJ). 
For the Fr atom, we use a large uncontracted [9s9p9<i] Gaussian basis set yielding 162 
molecular orbitals (table EJ). Molecular calculations for RbFr and CsFr then combine 
the Fr basis with one of the above basis sets for Rb and Cs. 

Core-valence correlation is treated by adding a phenomenological core polarization 
operator Wqpp to the Hamiltonian of the valence electrons [13J in the field of the M + 
and M' + ionic cores, and : 

W CP p = --a d f u - ^d Jm ( 5 ) 



where aj is the static dipole polarizability of the M + ion, and f\[ (resp. f' h 



M) 



is the total electric field seen by M + (resp. M' + ) due to the two electrons i and the 
partner ion M' + (resp. M + ): 



V^ ^iM r^, \ , Zm'Rm'M 



fu = 2^ — —9{riM) + —^ (6) 



=1,2 ' iM P^M'M 
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Atom 


Basis 


£ 


Exponents 


Cs 


[9s6p6d4:f] 


s 

P 
d 
f 


2.35,1-492561,0.824992,0.5,0.234682,0.032072,0.013962,0.005750,0.0025 
0.22,0.128,0.040097, 0.03,0.014261,0.004850 
0.29,0.12,0.096036,0.026807,0.009551,0.004 
0.2,0.1,0.05,0.005 


Fr 


[9s9p9d] 


s 

p 
d 


1.8,0.480468,0.369521,0.2,0.11230,0.053409,0.018240,0.006,0.002 

0.22,0.120,0.0655,0.03,0.0162,0.008,0.00443,0.002,0.001 

1.3,0.6,0.30,0.196894,0.067471,0.027948,0.010712,0.00300,0.001 



Table 4. Exponents of the Gaussian functions introduced in the [9s6p6dAf] basis for 
Cesium (basis "B"'), and the [9s9p9d] basis for Francium. 



Atom 


Series 


«r 


Ps 


Pp 


Pd 


Fr 


A 


20.38 [62J 


3.16372 


3.045 


3.1343 




B 


23.2[53 


3.343 


3.435 


3.292 


Rb 


A 


9.245 


2.5213 


2.279 


2.511 




B 


9.245 


2.5538 


2.3498 


2.5098 


Cs 


A 
B' 


15.117 
16.33 


2.6915 
2.0081 


1.8505 
2.6865 


2.807 
2.83518 



Table 5. Dipole polarizabilities a^ 1 and cutoff parameters pg introduced in the 
effective core polarization potential (CPP) for M + =Fr, Rb, Cs. Values for Rb (basis 
A and B) and Cs (basis A) are recalled from paper I, while new values for the Cs B' 
basis are reported. 



Experimental values for a^ are available for Cs and Rb as reported in paper I. 
We introduce a dependence on the electronic orbital momentum P. in the effective core 
polarization potential [33] , through the cut-off functions O^ to prevent the CPP matrix 
elements from diverging at the origin. As in ref. [33], ©/ is taken as a step function 
vanishing below the ^-dependent cut-off radii pg. 

For each basis set and a fixed value of a^ 1 , the cut-off radii are adjusted to 
reproduce the experimental energies [HUES EH] of the lowest s, p and d levels of Cs and 
Rb (Table EJ). As no experimental determination for the Fr + polarizability is available, 

=20.38 a.u. [H23- Another estimate 



we use like in the previous MP approach a^ r 



a^ r =23.2 a.u. was obtained from an extrapolation of known core polarizabilities of 
the lighter alkali atoms [HS], and we will discuss in the next section the influence of 
this parameter. As for Rb and Cs, the cut-off radii are adjusted on the experimental 
atomic energies jSHl E£]- for each above value of af r (Table EJ). The sensitivity of the 
results with the cut-off radii values is such that the influence of the polarizability can 
be balanced to provide a good representation of energy levels with both sets A and B. 



Electronic properties of Francium diatomic compounds and prospects for cold molecule formationlO 



Level 


Exp 


with a% r+ =20.38 a.u. 


AE 


with a^ r =23.2 a.u. 


AE 


7s 


-32848.87 


-32848.87 


0.002 


-32848.87 


0.002 


7p 


-19487.07 


-19486.93 


0.143 


-19487.18 


-0.110 


6d 


-16471.07 


-16470.92 


0.147 


-16470.94 


0.128 


8s 


-13116.62 


-13161.02 


-44. 


-13136.97 


-20. 


8p 


-9372.35 


-9448.16 


-76. 


-9424.96 


-53. 


Id 


-8550.96 


-8394.77 


156. 


-8371.59 


179. 


9s 


-7177.85 


-7161.96 


16. 


-7151.28 


26. 


9p 


-5573.18 


-5619.67 


-46. 


-5605.27 


-32. 


8d 


-5221.39 


4981.02 


240 


4963.01 


258 



Table 6. Computed binding energies of atomic Francium atom in cm -1 with the 
two available Fr + polarizabilities, compared to experimental values |33l I4fi) . and 
corresponding energy differences AE. For the unobserved 9p level our reference is the 
spin-averaged value predicted by the model-potential approach described earlier in the 

text. 



3. Molecular calculations 



As in paper I, we used an automatized procedure based on the CIPSI package [71] . The 
molecular orbitals are determined by restricted single electron calculations including 
the CPP [HI], namely corresponding to [XY] + systems, providing the potential curves 
for the relevant molecular cations. A full valence configuration interaction (CI) is 
then performed for each involved molecular symmetry, providing potential curves 
and permanent and transition dipole moments. In the case of neutral pairs, two- 
electron full CI is achieved in the framework of a single-component scheme. Clearly, 
further investigations using two- and four component schemes or larger cores would be 
interesting, although not easily tractable even at present, especially in the case of excited 
states, if one wants to include excitations corresponding to core polarization effects. Such 
calculations would anyway imply the use of two- or four-component pseudopotentials, 
which was not achieved here. 

In addition to the explicit above calculations, it is also important to estimate the 
magnitude of the deviation from the Coulomb repulsion of core-core interactions, which 
could be significant for the heavier and larger alkali atoms. The most significant effect 
is the repulsive overlap of the cores, which can be either extrapolated from exponential 
formulas valid for the lighter species [14], or derived from a calculation of [XY] ++ at 
the Hartree-Fock level, for instance within the frozen core approximation [T7[ IT£| 172*] 
which is adapted here as the core polarization is accounted for through an effective 
potential. The main influence of this term is to make the short-range repulsive part 
of the potentials steeper, while it dies out in the region of the minimum of the ground 
state. According to ref.[Z2], this term contributes for about 15 cm" 1 at the equilibrium 
distance of the ground state of Na 2 , K 2 , Rb 2 and Cs 2 , and slightly more (ps40 cm -1 ) for 
RbCs. A comparable estimation is probably appropriate for the Francium compounds 
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studied here. 

The contribution of the core-core dispersion energy can be estimated via the London 
formula (see for instance ref.|14j): 

v dis P (T}\ _ 6a d a d -frj f7 ( 7 \ 

cc { ' 2R* Ef + +Ej + {) 

where Ef and Ej are the ionisation energies of the X + and Y + ions respectively. 
At the equilibrium distance of the ground state, this term deepens the potential by 
about 27, 16 and 18 cm -1 for Fr2, RbFr and CsFr respectively. It basically compensates 
the core-core term above. 

3.1. Potential curves of Fr\ , RbFr + , CsFr + 

We display in figure El the potentiel curves of the ground state (1) 2 £+ and the lowest 
(1) 2 E+ state of Fr^~. We see in figure^ that the choice of the large Fr + polarizability 
af r =23.2 a.u. induces a potential well deeper by about 16 cm -1 than with 
a^ r =20.38 a.u. . Its effect is then rather small, possibly within the usual inaccuracy 
which was previously observed for this type of quantum chemistry calculations. 

Since no experimental result is presently available, we qualitatively compare our 
results for the Fr^" ground state with those of Rb^~ and Cs^~ in the main panel of Figure 
H3 and in Table EJ The well depth of the ground state indeed decreases for increasing 
mass of the atom, as it can be observed over the whole sequence of homonuclear alkali 
ions. In contrast the Frj" equilibrium distance is found slightly smaller the Cs^ one, and 
larger than that of Rb^". Note that a similar trend was also found by Lim et al jH] , but the 
large dispersion of the dissociation energy found by these authors does not allow an easy 
comparison with our results. These results illustrate that the properties of the Fr atom 
are not reducible to extrapolations from the lightest ones due to relativistic effects. This 
is fingerprinted by the size evolution of the lowest ns valence orbital, the 7s orbital of 
Francium being more contracted than the 6s one of cesium, whenever relativistic effects 
are taken into account. Indeed, Dirac-Fock calculations yields averaged radii values of 
5.91 ao for the 7s orbital in Francium and 6.08 ao for the 6s orbital in cesium, while 
the non-relativistic HF calculation reverses this hierarchy yielding 6.25 ao for Francium 
7s, and 5.85 ao for Cesium 6s. Such orbital contractions were illustrated on atoms by 
Desclaux [ZQj. In the case of alkali's, these contractions are further correlated with an 
increase of the atomic ionization potentials as discussed later on in the paper. This 
behaviour illustrates on francium compounds the influence of relativistic effects on the 
chemical bonding whenever heavy atoms are involved. 

A similar pattern is observed for the (1) 2 £+ state (jBb). It is interesting to note 
that such a shallow long-range potential well is predicted for all ionic alkali dimers 
including heteronuclear ones [731 EH US] i an d could be used as an alternative route for 
the detection through resonant ionization of ultracold molecules (121 EH! • 

We performed the same analysis for the two lowest states of RbFr + and CsFr + 
compared with those of RbCs + (figure EKnd Table EJ. Among the three atoms, Rb 
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Figure 3. Main panel: potential curves of the ground state 1 2 S+ (full lines) and the 
lowest 1 2 S+ state (dashed lines) of Fr^~ (black lines), compared with those of RbJ 
(red lines with open circles) and Cs^~ (blue lines with closed circles). All curves are 
referenced to the same dissociation energy taken as the origin, (a) Comparison of 
calculations with a^ r =20.38 a.u. |^2] (full line) and with a^ r =23.2 a.u. |53~] 
(dot-dashed line), (b) Blow up of the long-range well predicted in the 1 2 £+. The main 
and (b) panels display calculations with a^ r =20.38 a.u., and with basis B and B' 
for Rb^ and Cs;j~ repectively. 

Figure 4. Potential curves of the ground state (1) 2 £ + (panel a) and of the (2) 2 E + 
state (panel b) of RbFr + (black lines) and CsFr + (full red lines with open circles), 
compared with those of RbCs + (blue lines with closed circles). The energy origin is 
taken in both cases at the dissociation limits of the potentials. Calculations reported 
in panels a and b are performed with a^ r =20.38 a.u. [H2J, and basis B and B' 
for Rb and Cs respectively, (c) Ground state of the RbFr + ion: calculations with 
a^ r =20.38 a.u. and basis A for Rb (full line), with a^ r =23.2 a.u. jSSJ and basis 
A (Rb) (black dot-dashed line), and with a^ r =20.38 a.u. and basis B (Rb) (full line 
with circles), (d) Same as (c) for CsFr + ion. 



has the largest ionization potential and Cs the lowest one. The ground state (1) 2 £ + of 
RbFr+ (resp. CsFr + ) is correlated to Rb(5s)+Fr+ (resp. Fr(7s)+Cs + ) and the (2) 2 £+ to 
Fr(7s)+Rb + (resp. Cs(6s)+Fr + ). There is no clear trend for the hierarchy among well 
depths of the three systems. An analysis of the ground state of the whole ensemble of 
heteronuclear alkali molecular ionic species [71] shows that the NaLi + has the deepest 
well (more than 8000 cm -1 ) with the smallest equilibrium distance (around 6.5ao). The 
RbCs + , KCs + and KRb + ions have comparable well depth (in the 5000-6000 cm -1 
range) and equilibrium distances (around 8.5-9ao), while the remaining species exhibit 
well depth ranging between 3500 and 5000 cm -1 , with comprable equilibrium distance 
around 8-9ao- The RbFr + and CsFr + ions fit into the second category. As already 
found for RbCs + in ref.jZH], a shallow long-range well is present in the (2) 2 £ + potential 
curve of RbFr + and CsFr + ions (figure 0]d) similar to the well of (1) 2 S+ potential of 
homonuclear dimers. The difference induced by the different choices of parameters 
within these compounds is illustrated in figure Efc,d. As in Fr^", the larger Francium 
polarizability yields a deeper potential well. Increasing the size of the Rb basis has 
the same effect, as the RbFr + well depth is increased by about 11 cm -1 . In contrast, 
increasing in the same way the size of the Cs basis slightly changes the CsFr + equilibrium 
distance, but not its depth. 

3.2. Potential curves and permanent dipole moments for Fr2, RbFr, CsFr lowest states 

As no experimental data are available, the neutral systems involving Fr are analyzed 
along the same lines than the molecular ions in the previous section. Figure El shows 
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(i) 2 £+ 


Basis 


D e (cm- 1 ) 


R e (a ) 


Fr+ 


(a) 
(b) 


5537 

5552 


9.61 
9.65 


Rb^ 


A 
B 


6138 
6208 


9.08 
9.05 


Cs+ 


A 
B' 


5950 
5977 


9.75 
9.85 


((i) 2 s+ 


Basis 


D e (cm- 1 ) 


Re (Oo) 


Frt 


(a) 
(b) 


72.3 
72.6 


24.0 
24.0 


Rb+ 


A 
B 


81.7 
82.1 


23.0 
23.0 


Cs+ 


A 
B' 


82.7 
84.6 


24.3 
24.3 


(1) 2 £+ 


Basis 


D e (cm- 1 ) 


Re (Oo) 


RbFr+ 


A 
B 


5439 

5450 


9.3 
9.3 


CsFr+ 


A 
B' 


5054 
5052 


9.68 
9.72 


RbCs+ 


A,A 

B,B' 


5053 
5090 


9.40 
9.41 


((2) 2 S+ 


Basis 


D e (cm- 1 ) 


Re (Go) 


RbFr+ 


A 
B 


205.9 
206.1 


19.65 
19.65 


CsFr+ 


A 
B' 


291.3 
295.7 


19.07 
19.06 


RbCs+ 


A,A 

B,B' 


381.6 
389.5 


17.83 
17.84 



Table 7. Well depth D e and equilibrium distance R e (a =0. 0529177 nm) of 
the potential wells obtained in the present work for the two lowest 2 S states of 
Fr ionic compounds. Results for Rb^, CsJ and RbCs + are also displayed for 
comparison, (a) and (b) results for Frj~ correspond to a^ r =20.38 a.u. [H21 and 
with af r =23.2 a.u. [HSj respectively. Results for heteronuclear systems are obtained 
with a^ r =20.38 a.u.. Note that results with basis B' for Cs^ and (B,B') basis results 
for RbCs + are new, while those RbJ were already given in ref.|llj. 



the lowest Fr 2 potential curves dissociating into 7s + 7s, compared with those of Rb 2 
and Cs2- The influence of the value of the Fr + polar izability is clearly larger than in the 
molecular ion, as two electrons now contribute to the polarization of the cores (Figure 
EjD,d): the depth obtained with a^ r =23.2 a.u. [SH] is larger by about 80 cm -1 , and 
by about 12 cm" 1 for the (l) 1 ^^ ground state and for the (1) 3 S^ state respectively. 
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Figure 5. Potential curves of the ground state (1) 1 E+ (panel a) and the lowest (1) 3 E+ 
state (panel b) of Fr2 (full line with circles) calculated with a^ r =20.38 a.u. |62| . 
compared with those of Rb2 with basis B (blue dot-dashed line) and CS2 with basis B' 
(red dashed line). The energy origin is taken at the dissociation limits of the potentials, 
(b) Ground state of the Fr2 dimer obtained with with a^ r =20.38 a.u. (full line), 
and with a^ r =23.2 a.u. [Sj (dashed line), (d) Same as (b) for the (1) 3 E+ state. 

Figure 6. Potential curves of the ground state (1) 1 E + (panel a) and the lowest (1) 3 E + 
state (panel b) of RbFr (full black line) and CsFr (dashed red line), compared with 
those of RbCs (blue dots). The energy origin is taken at the dissociation limit of the 
potentials. We used a^ r =20.38 a.u. and basis B and B' for Rb abd Cs respectively. 

As for molecular ions the the well depth increases with increasing mass, while and the 
equilibrium distance increases along the series Rb2, Fr 2 , Cs2, again in agreement with 
ref.p. 

The lowest potential curves of RbFr and CsFr are displayed in Figure El compared 
with RbCs ones. In contrast with the related molecular ions, the RbCs ground state 
and lowest triplet state are found deeper than those of the Fr compounds. This suggests 
that the overall energy associated to the core polarization strongly depends not only on 
the polarizability of the individual atoms, but also on the difference of their ionization 
potentials, in a somewhat irregular way along the series of mixed alkali dimers. The 
characteristics of these wells are reported in TablelHl including those for Rb2, Cs2, RbFr, 
CsFr, and RbCs. As in the previous section, Rb2 results of type B were already shown 
in paper I, while those for Cs2 and RbCs are new. 

The permanent dipole moment of RbFr and CsFr are given in Figure [7| for the 
ground state and the lowest triplet state, respectively and compared with the results 
from paper I for RbCs. Note that results are not significantly different when using bais 
A or B and B' for Rb and Cs.The sign of charge transfer in the ground state of an XY 
diatomic molecule can be roughly determined by comparing the asymptotic energy costs 
A(XY) = IP(X) - EA(Y) versus A(YX) = IP(Y) - EA(X) describing ionization 
(IP stands for ionization potential) combined with electron attachement (EA stands 
for electron affinity) at dissociation. For alkali atoms from Li to Cs, the heavier the 
alkali atom, the smaller the ionization potential (43487, 41449, 35009, 33691, and 31406 
cm -1 ) and the larger the electron affinity (experimental values are 4984, 4403, 4043, 
3020, and 3803 cm -1 ), for Li, Na, K, Rb and Cs, respectively [49J. It is easy to show 
that in heteronuclear alkali diatomics including atoms up to Cs, the electron transfer 
from the heavy atom towards the light one is expected to be more favorable, due to the 
regular decrease of the IP's with increasing atom size which is not compensated by the 
corresponding variation of electron affinities. This is indeed confirmed for instance in 
our results from paper I, where we chose the convention of a negative sign of the ground 
state dipole moment in XY molecules corresponding to a X~Y + configuration where 
the X atom is lighter than the Y atom [T2*] . This model is no longer valid for Fr. Indeed, 
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(l) 1 ^ 


Basis 


D e (cm- 1 ) 


R e (a ) 


Fr 2 


(a) 
(b) 


3498 
3576 


8.45 
8.45 


Rb 2 


A 
B 


3861 
4039 


7.9 
7.9 


Cs 2 


A 
B' 


3703 

3787 


8.65 
8.69 


(i) 3 s+ 


Basis 


D e (cm- 1 ) 


Re (»o) 


Fr 2 


(a) 
(b) 


188.5 
200.8 


12.5 
12.5 


Rb 2 


A 
B 


233.6 

258.9 


11.6 
11.6 


Cs 2 


A 
B' 


255.0 
286.8 


12.0 
12.0 


(1)*E+ 


Basis 


D e (cm- 1 ) 


Re (Oo) 


RbFr 


A 
B 


3654 
3690 


8.2 
8.2 


CsFr 


A 
B' 


3553 
3576 


8.55 
8.57 


RbCs 


A,A 

B,B' 


3655 
3921 


8.27 
8.29 


((1) 3 £+ 


Basis 


D e (cm- 1 ) 


Re (Go) 


RbFr 


A 
B 


207.2 
209.8 


12.0 
12.0 


CsFr 


A 
B' 


209.5 

217.9 


12.36 
12.33 


RbCs 


A,A 

B,B' 


253.7 

273.1 


11.70 
11.65 



Table 8. Well depth D e and equilibrium distance R e of the potential wells obtained 
in the present work for the ground state and the lowest triplet state of Fr neutral 
compounds. Results for Rb2, Cs2 and RbCs are also displayed for comparison, (a) and 
(b) results for Fr2 correspond to a^ r =20.38 a.u. [52] and with a^ r =23.2 a.u. |S5] 
respectively. Results for heteronuclear systems are obtained with a^ r =20.38 a.u.. . 
Note that results with basis B' for Cs2 and with (B,B') basis for RbCs are new, while 
those Rb2 were already given in paper I. 



the ionization potential IP (Fr) =32848 cm -1 is larger by 1442 cm -1 than the one of Cs, 
almost half-way between IP(Cs) and IP(Rb), whereas the electron affinity EA(Fr)=3920 
cm -1 (HI] (see also ref. [ZZj) exceeds EA(Cs) by only 117 cm -1 . This IP increase is 
clearly related to the decrease of the orbital radius mentionned above. The energy cost 
A(CsFr) is lower than A(FrCs) by 1559 cm -1 , while A(RbCs) is larger than A(CsRb) 
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Figure 7. Permanent dipole moment functions of the ground state (1) 1 S + (panel a) 
and the lowest (1) 3 E + state (panel b) of RbFr (black dashed line) and CsFr (red dot- 
dashed line, ompared with those of RbCs (blue full line). A negative dipole moment 
corresponds to a X~Y + configuration. We used a^ r =20.38 a.u. and basis A for Rb 
abd Cs. 

by about the same amount. The charge transfer in CsFr is then expected to be Cs + Fr~, 
which is confirmed by the positive sign of the dipole moment. The situation of RbFr is 
more ambiguous since the difference A(RbFr) — A(FrRb) = —57 cm -1 is quite small, 
predicting a Rb + Fr~ arrangement. But the dipole moment is found negative and of 
small magnitude negative, compatible with a very weak electron transfer to form a 
Rb~Fr + arrangement. This apparent disagreement may be due to an uncertainty on 
the calculated electron-affinity of Francium [ST], but it most probably demonstrates the 
role of the polarization forces occurring at finite distance, resulting from the subtle and 
self-consistent unscreening of the ion electric field by the one of the electrons. 

The equilibrium distance of the lowest triplet state is large, and the state behaves 
like a neural state due to spin-forbidden electron transfer. The dipole moment then 
remains very weak, but we remark that its variation for RbCs is opposite to the RbFr 
and RbCs ones. At shorter distances, the magnitudes of the dipole moments may be 
increased due to overlapping atomic distributions. 

3.3. Excited states of Fr 2 , RbFr, CsFr and prospects for photoassociation and cold 
molecule formation 

It is also interesting to look at the Fr 2 excited states dissociating to the 7s + 7p limits 
in order to verify if the Fr 2 system behaves similarly to the Cs2 and Rb 2 species. 
For instance the existence of long-range molecular wells [7H] has been decisive in the 
formation process of Cs 2 and Rb 2 molecules jZHUHOj- Excited potential curves of gerade 
and ungerade symmetry correlated to the 7s + 7p limit are reproduced in FigurelHl Their 
overall variation is similar to the one for all homonuclear alkali pairs, with for instance 
a crossing at short distances between the (l) 1 ^^ and (l) 3 iT M curves, which gives rise to 
an avoided crossing when spin-orbit interaction is introduced. 

The magnitude of the atomic spin-orbit coupling is reported in the figure, showing 
that it will be the dominant interaction down to quite short internuclear distances 
(R « 20a ). Assuming that a constant spin-orbit interaction is a reasonably good 
approximation for R > 20a (as it was the case in Cs 2 , see for instance ref. jHJJ), we 
can diagonalize the electronic Hamiltonian including spin-orbit (also given in ref. |81j), 
yielding the curves reported in Figure El In contrast with the Cs 2 |H2] and Rb 2 JEJJ 
molecules, the spin-orbit interaction is so strong in Francium that it dominates the 
electrostatic interaction over a range of internuclear distances such that it prevents 
the occurence of long-range potential wells in the 0~ and l u symmetries, resulting 
respectively from the interaction between the 3 Il g and 3 £j~ states, and between the 3 I1 U , 
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Figure 8. Excited potential curves of gerade (panel a) and ungerade (panel 
c) symmetry correlated to the Fr(7s)+Fr(7p) dissociation limit, obtained with 
a^ r =20.38 a.u... Panel b and d focus on the long-range part of these curves. 

Figure 9. Same as Figure 03 but including now the fine structure induced by the 
atomic spin-orbit interaction. 

1 I1 U , and 3 £+ states. For the same reason, the expected avoided crossing generated by 
the interaction between the (l) 1 !]^ and (1) 3 I1 U states is no more well localized so that 
the resulting 0+ curves are almost uncoupled. 

The photoassociation rate for cold Francium atoms can be estimated by comparison 
with the cesium one, assuming that similar experimental conditions could be achievable. 
Following the model developed in ref.jBHJ, we consider a PA transition close to an 
atomic D line. The photoassociation rate Rpa is then proportional to the product of 
characteristic constants A^/i -1 / 2 ^ r^ 1 , where Xd is the wavelength of the D transition, 
\i the reduced mass of the atom pair, C% the effective long-range coefficient of the PA 
state, and tjj the radiative lifetime of the D line. Then the ratio of PA rates is estimated 
according to: 

R PA (Fr) f\ D (Fr)\ 3 ( ^ Cs \ 1/2 [C[ r \ 2/3 r D (Cs) 



Rpa(Cs) \X d (Cs)J V/W \C§* r D (Fr) 



where the effective C3 coefficient for Francium has been roughly estimated 
from the atomic lifetime of the D2 line [12] ■ Equation |H1 results in a factor of about 0.7, 
meaning that PA is expected to be almost as efficient in Francium than in Cesium. On 
the other hand, the absence of double-well potential, or of a significant interaction 
between the two 0+ curves (labelled as resonant coupling in ref.|85j) will probably 
prevent the efficient formation of cold Francium dimers in deeply bound levels of the 
ground state or of the lowest triplet state (see refs. [23 EH] for a full description of 
these mechanisms), leaving molecules only in the uppermost bound vibrational levels. 
Moreover, the van der Waals coefficient of ground state Francium Cg =5226 a.u. is 
smaller than the one of Cesium C§ s =6851 a.u.jTUj. The corresponding dimer potential 
curves have a shorter range, which also does not favor the stabilization of long-range 
PA molecules by radiative decay towards these states. 

The excited potential curves correlated to the two lowest s + p dissociation limits 
in RbFr and CsFr are drawn in Figure E3 as obtained with a^ r =20.38 a.u., and 
with basis B and B' for Rb and Cs. In contrast with all other the mixed pairs, the 
smallest excitation energy is now for the lightest atom, which is another characteristic 
which illustrates the unusual character of the Francium atom compared to the rest of 
the alkali series. In RbFr, the spin-orbit is expected to strongly mix both bp + 7s and 
5s + 7p asymptotes (see the dashed and dotted boxes in Figure E3 which cannot be 
represented by the perturbative approach employed for Fr 2 , and multiple- well structures 
may probably be expected. The CsFr asymptotes are more spaced than in RbFr, so 
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Figure 10. Excited potential curves for (a) FrRb and (b) FrCs molecules correlated to 
the two lowest s+p dissociation limits, obtained with a^ r =20.38 a.u., and with basis 
B and B' for Rb and Cs respectively. The dashed and dotted area indicate respectively 
the atomic fine structure splitting for Rb and Fr (in panel a), and for Cs and Fr (in 
panel b). Full lines: singlet states. Dashed lines: triplet states. Black lines: £ states. 
Red lines with circles: n states. The energy origin is taken at the 5p + Is and 6p + 7s 
limit for RbFr and CsFr respectively. 

Figure 11. Transition dipolc moments from the ground state towards the first 
excited (a) 1 E + states and (b) 1 H in Francium compounds. Full black lines: X — » 
(l) 1 E+(7s + 7p) and X -> (l) 1 n„(7s + Ip) in Fr 2 . Full red lines with closed circles: 
X ->■ (l) 1 S+(5p + 7s) and X -v (l) 1 U ( ~5p+7s) in RbFr. Dashed red lines with closed 
circles: X -> (2) 1 S+(5s + Ip) and X -> (2) 1 n(5s + Ip) in RbFr. Full blue lines with 
open circles: X -* (l) 1 £+(6p+7s) and X -» (l) 1 n(6p + 7s) in CsFr. Dashed blue lines 
with open circles: X -> (2) 1 S+(6s + 7p) and X -» (2) 1 n(6s + 7p) in CsFr.Calculations 
are performed with a^ r =20.38 a.u., and with basis A for both Rb and Cs. 

that the perturbative approach may be valid at large distances, but strong mixing will 
probably occur at short interatomic distances. 

The van der Waals coefficients for these systems, C§ bFr =4946 a.u. and 
qCsFt = 5ggg au have the same magnitude than the RbCs one CQ bCs =5663 a.u. [333], 
so that the stabilization of PA molecules by radiative decay should be as efficient as in 
RbCs |35| . 

Finally we present in Figure E] the transition dipole moments between the ground 
state and the first excited 1 S + states and 1 I1 in Francium compounds. Here also results 
are equivalent when using basis A or basis B and B' for Rb and Cs respectively. In 
the upper panel, one sees that the RbFr and CsFr dipole moments for the transition 
towards the lowest £ excited state are superimposed onto the Fr 2 function between 7 
and 12ao, suggesting that in this range, the electrons are closer to the Francium cores 
confirming the discussion about permanent dipole moments above. The same pattern 
is observable in panel b as well, but over a shorter range (8 to 10 do). The presence of 
a strong avoided crossing between the 1 I1 states in RbFr and CsFr is visible here as the 
two transition dipole moment functions cross each other at the position of the avoided 
crossing. 

4. Conclusion 

In this work we investigated the possibility to create cold Fr 2 , RbFr, and CsFr molecules 
through photoassociation of cold atoms. Indeed, the creation of dense samples of cold 
Francium atoms may become available experimentally, as well as the mixing with cold 
Rubidium or Cesium atoms. We predict that under comparable experimental conditions, 
photoassociation of cold Francium atoms should be as efficient as with cold cesium 
atoms. In contrast the formation of ultracold Fr 2 molecules will probably be weaker as 
the formation processes known in Cs 2 , relying on the presence of long-range potential 
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wells or resonant coupling, are not expected in Fr 2 . However, the formation of ultracold 
RbFr or CsFr molecules is probably as efficient as in RbCs, which could be worthwhile 
to try experimentally as dense cold Rb and Cs samples are widely available. 

To extract these results, we performed the first calculations of the electronic 
structure of the Fr2, RbFr, and CsFr molecules, yielding potential curves, permanent and 
transition dipole moments. A new scalar (ie spin-orbit averaged) pseudopotential has 
been determined to represent the interaction of the valence electron with the large Fr + 
core complemented with core-polarization operators. Through a comparison with the 
closest heavy alkali dimers Rb 2 , Cs2 and RbCs, we found that our results yield a realistic 
representation of the Fr2, FrCs and FrRb molecules. It would be interesting to check the 
present results with two-component pseudopotential configuration interaction schemes 
|87| . but they only exist in the framework of small core pseudopotential (ie 10 valence 
electrons or more in he case of Francium. This would allow in particular to estimate 
errors of the present calculation which ignore spin-orbit induced polarization effects. 
However, we provided a description of the molecular fine structure via an empirical 
atom-in molecules schemes, which suggests that the behaviour of excited Fr 2 molecular 
states involved in photoassociation and cold molecule formation is different from Rb 2 
and Cs2 cases. Here again, it would be interesting to compare with more refined ab 
initio spin-orbit calculations which are in progress. An interesting conclusion of the 
present work is that, while the nature and the quantitative behaviour of the ground state 
and metastable triplet state potential curves of the francium compounds are generally 
comparable with other previously studied alkali diatomics, the ground state polarity of 
hetero-atomic molecules cannot be extrapolated from the other alkali pairs, due to the 
influence of relativistic effect which affect atomic properties, especially the ionization 
potential(correlated with the outer s orbital radius). We hope that the present work will 
stimulate experimental spectroscopic studies on those short-lived species, which would 
also help to discuss further the accuracy of the present results. 
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